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We describe criteria for implementation of quantum computation in qudits. A qudit is a d-dimensional system 
whose Hilbert space is spanned by states |0), 1), . . ., \d— 1). An important earlier work of Mathukrishnan and 
Stroud [ I ] describes how to exactly simulate an arbitrary unitary on multiple qudits using a 2rf — 1 parameter 
family of single qudit and two qudit gates. Their technique is based on the spectral decomposition of unitaries. 
Here we generalize this argument to show that exact universality follows given a discrete set of single qudit 
Hamiltonians and one two-qudit Hamiltonian. The technique is related to the QR-matrix decomposition of 
numerical linear algebra. We consider a generic physical system in which the single qudit Hamiltonians are a 
small collection of H* k = h£l(\k){j\ + \j){k\) and = h£l(i\k){j\ — i\j){k\). A coupling graph results taking 
nodes 0, . . ., d — 1 and edges j <-> k iff H x ^ are allowed Hamiltonians. One qudit exact universality follows 
iff this graph is connected, and complete universality results if the two-qudit Hamiltonian H = —ftQ.\d — 1 , d — 
l){d — l,d — 1| is also allowed. We discuss implementation in the eight dimensional ground electronic states of 
87 Rb and construct an optimal gate sequence using Raman laser pulses. 

PACS numbers: 03.67.-a, 03.67.Lx 



I. INTRODUCTION 

An important theoretic construct used in the field of quan- 
tum information is the qubit. Its utility follows from the sim- 
ple but significant recognition that all two dimensional sub- 
spaces, regardless of the underlying physical system, can be 
regarded as informationally equivalent. This has made it pos- 
sible to discuss quantum computation in terms of single qubit 
and two qubit gates without the need to analyze the specific 
interactions that realize operations within a physical system 
or between subsystems. An important issue in this regard is 
that a necessary condition for efficient quantum computation 
is the existence of an underlying tensor product structure on 
the Hilbert space H . If all computation were performed on a 
single d — dim{fK) level system then some physical resource 
such as space or energy would grow with the dimension of the 
system [2]. In contrast, the analogous resources grow poly- 
logarithmically with the dimension when the system is com- 
posed of many subsystems. By this argument, a computation 
performed on qubits (d = 2) is in some sense the most effi- 
cient foliation of Hilbert space. 

Nevertheless, there are compelling reasons to consider 
computation on qudits with d > 2. First, most physical im- 
plementations encode qubits in a subspace of a larger Hilbert 
space. Using higher dimensional subspaces already endowed 
in these systems may be more efficient in terms of the num- 
ber of interacting gates needed for an algorithm that acts on a 
Hilbert space of fixed dimension. This is critical for error con- 
trol because interactions between qudits tend to open channels 
for interactions with the decohering environment. By contrast, 
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in many physical systems, single qudit control is a well devel- 
oped technology that can be done with high precision. Sec- 
ond, there is some evidence that the error thresholds for fault 
tolerant computation improve when the encoding is done with 
qudits where d > 2 and prime [3]. 

Previous work by Brylinski and Brylinski proves the nec- 
essary and sufficient criteria for exact qudit universality [4]. 
Exact universality means that any unitary and, by unitary ex- 
tension to a larger Hilbert space, any quantum process, can 
be simulated with zero error. The authors show that arbitrary 
single qudit gates complemented by one entangling two qudit 
gate is needed. Their method is not constructive. Muthukrish- 
nan and Stroud [1] give a constructive procedure for an exact 
simulation of an arbitrary unitary on n qudits using single qu- 
dit and two qudit gates. Their approach uses the spectral de- 
composition of unitaries and involves a gate library consisting 
of a family of continuous parameter gates. Here we describe 
an approach that uses the QR decompositions on unitaries to 
achieve exact universal computation on qudits. This construc- 
tion has the advantage that the single qudit gates are generated 
by a fixed set of Hamiltonians that couple pairs of states in the 
single qudit logical basis. The gates perform rotations, param- 
eterized by one angle, about orthogonal axes within the asso- 
ciated two dimensional subspace. Additionally, our decom- 
position requires only one fixed two qudit gate, the controlled 
increment gate (CINC) gate. This gate can be simulated by 
at most d — 1 instances of a two-qudit Hamiltonian that 
generates a phase on a single product state of two qudits. Such 
interactions can be engineered in many atom optical systems. 

In this paper, the general results are developed with close 
contact to the example of computation in the d = 8 qudit en- 
coded in the ground hyperfine states of 87 Rb. In Sec. II we 
describe the construction of single qudit unitaries using the 
QR decomposition. We introduce a coupling graph to describe 
how states are connected to each other by physical Hamilto- 
nians. The set of rotation planes may be incomplete, i.e. each 
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state may not be connected to every other state. However, 
provided the graph is connected, an efficient decomposition 
can be found. Multiqudit computation is addressed in Sec. 
III. It is shown that a single two-qudit gate when combined 
with single qudit gates suffices to generate arbitrary two qu- 
dit unitaries and hence completes the requirements for exact 
universality. In the appendix we show how to convert between 
the gate library introduced here and the family of gates used in 
Ref. [1]. Finally, we conclude with a summary of the results 
in Sec. IV. 



H. ONE-QUDIT UNITARIES 

We pick a fixed gate library for single qudit operations in- 
volving rotations about non orthogonal axes of two dimen- 
sional subspaces. Within each subspace Hn, the gates are 
generated by the two Hamiltonians: 



H% = h£l(\j){k\ + \k)(j\), 



H 



■ h£i(-i\j)(k\ + i\k)(j\). 



(1) 

For convenience of notation, we assume the strength of each 
coupling is equal to £2 and leave the time each Hamiltonian 
is applied as a free parameter. Any unitary in the two dimen- 
sional subspace can be written 



Uj fi(y,$,B) = exp 
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for the appropriate t,t',t" using the XYX Euler angle decom- 
position [5]. In some cases, the two Hamiltonians in Eq. 
1 can be turned on simultaneously. By adjusting the rela- 
tive strengths of the couplings, one can then realize any ro- 
tation about an axis on the equator of the Bloch sphere in the 
two dimensional subspace. For brevity, we write Ujk{y,§) = 
e -i(co S («^+sin(«ffJ)Y/(/iQ) ; where ^ . s understood that tf the 

couplings Hj k , H y - k cannot be turned on together then Ujt (y, <j)) 
requires three elementary gates. 

Realization of an arbitrary unitary evolution v € U(d n ) re- 
quires two steps. The first corresponds to a QR decomposition 
[6] of the matrix v. 

• Using the allowed set of Hamiltonians, we may re- 
alize matrices of Givens rotations physically. These 
are described in the next paragraph. Generically, the 
QR decomposition writes an invertible G = UT, where 
U = G\G2 ■ ■ ■ G( is a product of Givens rotations and 
hence unitary and T is upper triangular. Note that if 
G = V is unitary, then so likewise is T = U^G, whence 
T is in this case a diagonal matrix which applies relative 
phases to computational basis states. 

• Using techniques for realizing diagonal computations 
[7], a sequence of Hamiltonians realizing T is con- 
structed. 



We next illustrate the idea of a Givens rotation by way of 
example, retaining V as above. We may choose U a Givens 
rotation so as to zero the matrix element (UV)d-io (where 
the indices run 0, 1 , • • ■ , d — 1.) Specifically, suppose that 
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cos(y) — /e i(|) sin(y) 
-/e~ i(|) sin(y) cos(y) J 



Here, we choose the angles y, (]) as follows: 

tany = |v rf _i j0 /v rf _ 2j o| 

§ = 7c/2 + arg(v,i_2,o)-arg(vrf-io) 



(3) 



(4) 



where v m<n are the entrees of the unitary V. Then letting 
denote a changed entry, we obtain: 
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In the next step, one chooses a unitary U c i-\,d-3 to zero the 
matrix element (Ud-i d-2V)d-2fi- Continuing carefully in this 
way allows one to complete the QR decomposition described 
above by introducing a zero into every entry of the resulting 
unitary below the diagonal. 



A. Example: One-qudit Unitaries in 87 Rb 

We begin by describing explicitly the implications of our 
constructions for an example which is related to but not cov- 
ered explicitly by earlier work of Mathukrishnan and Stroud 
[ 1 ] . Specifically, we describe the coupling graph alluded to in 
the introduction in this case before defining it in general. 

Thus, consider the atomic species 87 Rb per Fig. 1. There 
are two ground state hyperfine manifolds with total spin = 
1 and Fi- = 2 split in energy by the hyperfine interaction E%t. 
Each manifold consists of 2F + 1 degenerate magnetic sub- 
levels Mp for a total of eight distinguishable states. We des- 
ignate this d = 8 system a quoctet. The degeneracy can be 
lifted by applying a longitudinal magnetic field B z . For small 
fields, the resultant Zeeman interaction is linear in the mag- 
netic quantum number: Hb = gFB z Mf, where the Lande g 
factors satisfy gp, = — gp^ [8]. 

There are several ways to couple the magnetic sublevels in- 
cluding the use of microwave pulses and Raman lasers. These 
techniques are usually distinguished by the strength of the 
coupling with respect to the hyperfine interaction. We con- 
sider coupling that is weak relative to Ef,f using a pair of laser 
beams on Raman resonance between two sublevels at a time. 
The effective atom-laser Hamiltonian Hal in the subspace Hi^ 



H A Ljk = cos($)H% + sm(<b)Hl 



(6) 
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Mp = -2 -1 1 2 

FIG. 1: A single d = 8 qudit (quoctet) encoded in the ground state 
hyperfine levels of 87 Rb. A pair of lasers can couple states in differ- 
ent hyperfine manifolds according to the selection rule AMf = 0, ±1. 
Projective measurements of population in state |7) are made by ob- 
serving resonant fluorescence on a cycling transition to the excited 
state. Any pair of states can be coupled by swapping neighbors to- 
gether pairwise and similarly any state can be measured by swapping 
to |7). 



where Q. = \Q.\Q.2\/A is the product of the individual laser 
Rabi frequencies divided by the detuning A from the ex- 
cited state, and 4> = <i>i — 4*2 is the relative phase of the two 
beams. In order to selectively couple two states it is necessary 
that their energy difference be unique. In the linear Zeeman 
regime, this can only be accommodated when the two lev- 
els reside in different hyperfine manifolds. Additionally, it 
will be important to minimize spontaneous emission during 
the pulse sequence by choosing a large detuning A of each 
laser from the excited states. The allowed couplings are con- 
strained by angular momentum selection rules which dictate 
the change in magnetic spin quantum number during a sin- 
gle pulse sequence. For detuning A much greater than the ex- 
cited state hyperfine structure, but less than than the fine struc- 
ture splitting, the angular momentum selection rules dictate 
AMf = 0, ±1 and |F|_,0) |-Ft,0). Using two-laser pulses of 
the appropriate frequency and polarization, the states \F\,Mp) 
and \Fj,Mf +AMf) where AMp = 0,±1 can then be coupled 
together. This is shown schematically in Fig. 1 where states 
|2) and |5) are coupled by a o + — Jt polarized laser pair. 

At this point we pause to comment on the resources neces- 
sary for single quoctet compution using Raman pulses. Tran- 
sitions realizing AM p = 0, ± 1 can be achieved by choosing the 
correct polarizations for the lasers with respect to a quantiza- 
tion axis defined by the magnetic field direction. For a fixed 
Zeeman splitting, it will be necessary to have lasers tuned to 
Raman resonance for all the allowed couplings. This may be 
achievable using a fixed source laser source that is frequency 
modulated appropriately. Another recourse is to change the 
magnetic field strength for each pairwise state coupling so that 
only one laser pair of fixed frequency is necessary. The phase 
shifts accumulated on the basis states during the change in 
Zeeman interaction can be accounted for in the gate sequence. 
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FIG. 2: This is the coupling graph for the coupled hyperfine states 
of 87 Rb (cf. Fig. 1.) As it is connected, the collection of atom-laser 
couplings allows for universal one-quoctet computation. 

We wish to show that the above set of atom-laser Hamil- 
tonians suffices to construct an arbitrary unitary evolution of 
the one-quoctet phase space H\ = C|0) © • • • © C|7). Take 
V € t/(8) as the target one-quoctet evolution, where t/(8) is 
the symmetry group of the inner-product on the Hilbert space 
(i.e. VV^ — The goal then is to decompose V into a se- 
quence of evolutions by these atom laser Hamiltonians: 

V = exp(iffj L *i !%)■■■ exp(iHi L t c /ti) (7) 

Additionally, we prefer efficient decompositions, i.e. we wish 
to use as few laser pulses (as small an I) as possible. This 
is sometimes not possible, depending on which states \ j), \k) 
are coupled by an H^l. In order to classify when the QR step 
is possible, we introduce the notion of a coupling graph, by 
example. 

87 Rb coupling graph: The 87 Rb coupling graph has vertices 
labelled by 0, 1 , . . . , 7. In addition, consulting Fig. 1, we also 
allow in the following edges, corresponding to the atom-laser 
coupled hyperfine states. 

{(0,5), (0,6), (0,7), (1,4), (1,6), (2,3), (2,4), (2,5)} (8) 

In particular, the edges encode the selection rule for the hy- 
perfine states. The graph is reproduced in Fig. 2. We note for 
future use that it is connected. Provided the states \j), \k) are 
coupled, we may produce any determinant-one unitary evolu- 
tion of ttjk using Eq. 2. 

Now note that since the coupling graph is connected, we 
may in fact sequentially construct a Givens rotation on any 
ttfe. Indeed, even if \j) and \k) are not paired, there exists 
a sequence \j ) = \j), \ji), \j 2 ), . ■ ■ , \je) = \k) such that each 
consecutive pair admits atom-laser Hamiltonians. Moreover, 
taking (j) = 7t/2,0 = 7i/2in Equation 3 shows that we may use 
these pairings to swap states up to relative phase. Hence, since 
we may physically construct some sequence of Hamiltonians 
for any Givens rotation, we see that the first step of the QR 
decomposition is possible. 

This leaves open the question of efficiency. For example, 
one might hope that in a graph as highly connected as that for 
87 Rb few or no swaps might be required. This is indeed pos- 
sible as we now show. It is convenient to reorder the unitary 
in a logical basis labeled {7,0,6,5,3,2,4,1}. By successive 
Givens rotations, one may bring a unitary V to diagonal form 
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column by column where the sequence is chosen so as to not 
void zeroes created in earlier steps. Each of the columns can 
be reduced to a single unimodular entry on the diagonal by 
a sequence of Givens rotations Uj.k acting on the two dimen- 
sional subspace ttjk as follows [10]: 

• Column 7 reduction: t/^it/i^t^, 3^5,2^0,5^0.6^7 ,0 

• Column reduction: U4.1U2.4U2. 3U5, 2^0,5^0,6 

• Column 6 reduction: C/2,5C/2,3C/4,2C/i,4t^6.i 

• Column 5 reduction: 1/4^1/2,41/2,3^5,2 

• Column 3 reduction: t/4,1 1/2,41/3,2 

• Column 2 reduction: L/^i L/2,4 

• Column 4 reduction: i/4,1 

Note that in general, constructing Un requires 2d(j,k) — 1 
basic Hamiltonians, where d(j,k) is the distance between j 
and k in the graph corresponding to the pairing relation. For 
qudit computation in 87 Rb using Raman pulses, the graph 
is sufficiently connected so that the distance is never greater 
than one in the QR decomposition above. The are a total of 
8 x 7/2 = 28 gates in the reduction to diagonal form. Each 
gate Uj.k £ SU (2) has two parameters so this gives 56 parame- 
ters. An arbitrary u £ SU (d) requires d 2 — 1 parameters so the 
additional 7 parameters correspond to seven relative phases 
left on the diagonal. 

B. Relative Phases 

The goal of this section is to show that should the Hamilto- 
nian graph be connected and T = TfjZo \j) U\ ^ e a diago- 
nal element of U(d), then we may realize T with the allowed 
Hamiltonians H\, H y - k . In fact, we only need to construct T 

up to a global phase so we can simulate T' = e"™- 1 7\ We first 
note that although it is not explicitly an allowed Hamiltonian, 
we may for any (_/,£) -edge within the coupling graph simulate 
the effect of H z jk = HQ.(\j)(j\ - \k)(k\). Indeed, for any fixed 
angle y we have 

e- m J^ m) =C/ ; -, t (-7t/4,Jt/2)C/ M (Y,0)i7 M (7t/4 ) Jt/2). (9) 

The goal then is to find a sequence of z rotations that simulates 
T: 

f[«p(-iff? 4 n/h) = T (io) 
1=1 

Given that the coupling graph is connected, choose a subset 
S of d — 1 edges V- k = \j){j\ — that leave the graph 

connected. We can represent the elements of S as vectors 
in a d dimensional real vector space spanned by the or- 
thonormal vectors {ej}, i.e. X z j k = ej — e^- We then con- 
struct a [d — 1) x d matrix M out of the row vectors in S: 
M = {^0*0^1*1 > ■ --K-2k d _ 2 }- The appropriate timings tj in 



Eq. 10 necessary to simulate T' are given by solutions to 
the matrix equation MB = (j), where 9 = £l{to, . . .td-2} an d 
(j) = {(j)o, . . .<j) ( /-2}- The angle (j) c /_i = for the unitary T' by 
assumption. It is easily verified by Gaussian elimination that 
the dimension of the row space of M is d — 1, thus there is a 
unique solution to the vector 8. 

The result is that any diagonal unitary can be simulated up 
to a global phase using 3 x (d — 1) gates from the gate library. 
This sequence can be reduced by a factor of three if z rotations 
can be implemented directly without conjugation. Further, all 
the Hamiltonians H% are diagonal and hence commute, so z 
rotations that act on disjoint subspaces can be implemented in 
parallel using additional control resources. 

C. One-qudit universality for generic coupling graphs 

We found that for computation in the 87 Rb quoctet, a sin- 
gle qudit unitary could be brought to diagonal form using the 
fewest possible Givens rotations. This is not peculiar to that 
system but is in fact possible for any system with a connected 
coupling graph. 

Lemma II.l ([10]) Given a d-node coupling graph Q of al- 
lowed Givens rotations, then any U £ SU(d) can be brought 
to diagonal form using d(d — l)/2 allowed rotations if and 
only if Q is connected. 

Proof: Suppose Q is connected. Form any spanning tree 
for it, and renumber the nodes so that the path from node d 
(the root of the tree) to any node j passes through no node 
numbered lower than j; such a numbering can be constructed 
by successively deleting leaf nodes and numbering in order of 
deletion. (For 87 Rb, we formed the tree by breaking the edge 
between nodes 6 and 1 and used the logical basis ordering 
{7,0,6,5,3,2,4, 1}.) At the ;'th step (J = 1, . . . ,d- 1), create 
the tree Tj, rooted at node j, from the portion of the spanning 
tree defined by nodes j,...,d. (Note that Tj is connected due 
to the way we numbered the nodes.) Then, until only the root 
of Tj remains, choose a leaf k, use a rotation defined by its 
edge to eliminate element (k, j) of U, and delete node k from 
Ti . The result of applying these steps is an upper triangular 
matrix (and therefore, since U is unitary, a diagonal matrix) 
computed by using d(d—\)/2 allowed rotations. 

Suppose g is not connected and consider a matrix U £ 
SU(d) that has no zero elements. Choose an arbitrary node 
to call node 1 . Then we can at best eliminate all but one of 
the nonzeros in column 1 of the disconnected piece, but there 
is no allowed rotation that will eliminate the last nonzero. Re- 
peating the argument for each choice of node 1 , we conclude 
that we cannot reduce U to diagonal form using only allowed 
rotations. □ 



III. MULTI-QUDIT UNIVERSAILITY 

Suppose in addition to being allowed local Hamiltonians 
{H x - k } with a connected coupling graph, the physical system 
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also allows for a two-qudit phase Hamiltonian 

Hi m = -h£l\d-l,d-\){d-\,d-l\. (11) 

Using known but dispersed techniques [1, 11], we describe 
a bootstrap which allows for universal quantum computation. 
Note that due to the standard QR decomposition, it suffices 
to construct arbitrary two-qudit unitary evolutions [1,6]. In 
fact arbitrary one-qudit operations controlled on d— 1, i.e. 
Ai (v) = I d i_ d © v G U(d 2 ) for v G U(d), suffice. 

Before presenting the generic discussion, we describe a par- 
ticular example of a two-qubit operation which has seen heavy 
use [1]. First, we label as (©1) the self-map of TLjdTL which 
carries k i— » (k + 1 )mod d. Then the controlled-increment gate, 
abbreviated here as CINC, is defined by extending the follow- 
ing rule linearly: 

{ j.k J :\: V'!i ! • ™ 

The CINC gate is heavily used in the literature in building 
a generic ^-controlled computation A^(v) [1] as well as for 
constructing quantum error correction codes [12]. 

We may explicitly realize CINC from the Hamiltonian 
Hi,,, = — Q.\d — l,d — l)(d — l,d — 1| as follows. This dis- 
cussion uses the group theory notation that for ji,...,jt C 
{0,l,...,d — 1}, we write (7172 ■ ■■ji) for the cyclic permu- 
tation with ji i-> j z , h i3» • ' k-i jl> k ju and 
all other set elements fixed. The permutation will also be 
identified implicitly with the associated permutation matrix 

n Uih-Jt) e U ( d )- Hence ' S iven ( 01 )( 12 ) ■■■{d-2d-\) = 
©1, we see that CINC = Ai[(01)(12) ■■■{d - 2 d - 1)]. The 
construction of CINC then takes place in the following steps: 

• We may write exp(— ;'//,■„, Jt/(H2)) = Ai (I c /-2 © CJ Z ). 

• We next argue that we may construct Ai \{j j + 1 )]. In- 
deed, first note that using an appropriate single qudit 
permutation matrix Us ^, we may construct A] (Ij © a z © 

h-2-j) as 

M{I)®<f®Id-l-j) = /rf®f/y +2 ,rf-l(7t/2,0) 

Ai(7 rf _ 2 ©c^) 

/rf®f/y +2 ,rf-l(-7t/2,0). 

(13) 

Then 

Al[C/ ; + = /,/®f/;+l.;+2(n/4,7t/2) 

Ai(/j ffia J ffi/rf_2-;) (14) 

/d®t/;+l,;+2(-7t/4,7i:/2). 

• This leads to the realization of CINC in a max- 
imum of d — 1 controlled operations, given that 
CINC = Ai(01)Ai(12)---Ai(<*-2d-l). 

We finally consider the construction of an arbitrary Ai(v) 
for v G U(d). Again using standard Givens arguments, it 
suffices to construct Ai [Ij © w ®Id-j-i] f° r any w G U(2), 
det(w) = 1 . Indeed, using the block- wise permutation argu- 
ment above, A\(w®Id-i) suffices. Now recall ([1 1], Lemma 



5.1) that there exists for any w as above factor matrices a, b, 
and c such that w = a<3 x bo x c while I2 = abc. Hence 

Ai(wffi/ rf _ 2 ) = (^®a)Ai[(01)](7 rf (8i6)Ai[(01)](7 < ,®c) 

(15) 

This completes the bootstrap argument for exactly universail- 
ity, when the restricted one-qudit Hamiltonian set {Hjj*} is 
augmented by //„,, . 

We showed how the gate CINC can be constructed using 
the entangling interaction Hj„ t . In many situations, the inter- 
action between qudits will contain more than one term on the 
diagonal. For instance, the true Hamiltonian may be 

d-l 

H L= 52 h&mn\mn}{mn\. (16) 

mn=0 

In this case the interaction is entangling iff the following is 
true [4] 

Q. mn + Ci pq ^ Qmq + Cl pn + 2%k for any m,n,p,q, any k G Z. 

(17) 

When the interaction H\ nt is entangling, it is always possible 
to map it to //,„, using multiple applications of the coupling 
conjugated by single qudit gates. In practice, some multiqudit 
operations may be done more efficiently using H' jnt directly. 

There are several proposals for realizing diagonal coupling 
gates in real physical systems. For example, in trapped atoms 
possible coupling mechanisms include pairwise interactions 
via dipole-dipole interactions [13, 14], and controlled ground 
state-ground state collisions [15]. The later proposal has been 
realized recently between atoms trapped in an optical lattice 
[16]. These proposals were originally made with the goal 
of engineering two qubit controlled phase gates. As such, 
a naive adaptation to encoding over all magnetic hyperfine 
levels would fail due to off diagonal couplings between ba- 
sis states. However, it should be possible to modify one or 
more proposals to realize a differential shift on a single prod- 
uct state. For instance, in Ref. [ I 7] it was proposed to realize 
a quantum gate using the ground state-ground state collisional 
shift induced by shape resonance. Here one can tune a mag- 
netic field such that a single molecular state is on resonance 
with a bound motional state of an external trap for both atoms. 
Because the resonance is dependent on the internal states, a 
unique phase is accumulated on a single product state. Pro- 
vided the atoms are sufficiently separated, the other basis state 
pairs do not interact and a Hamiltonian of the form //„« is re- 
alized (up to local unitaries.) 



IV. CONCLUSIONS 

We have identified the criteria for exact quantum computa- 
tion in qudits. Our method is constructive and relies on the QR 
decomposition of unitaries on qudits using a gate library gen- 
erated by a fixed set of single qudit Hamiltonians and one two 
qudit entangling gate. Using the concept of a coupling graph 
we are able to show that universal computation is possible 
if the nodes (equivalently logical basis states) are connected. 
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Further we give a prescription for efficient single qudit com- 
putation by demanding that at each stage of the QR decom- 
position the graph remain connected. Using the gate library 
generated by the couplings in Eq. 1 the maximum number 
of gates is 3d(d + 1)/2 — 3. The technique for computation is 
exemplified with a quoctet using the Raman coupled magnetic 
sublevels of 87 Rb. It is shown that arbitary single quoctet com- 
putation is possible with at most k = A9 laser pulse sequences. 
This gate count is optimal and could be reduced to the mini- 
mum number k m i n = d(d — 1 )/2 = 28 only if one appends the 
diagonal generators H : - k to the library of couping Hamiltoni- 
ans. 

We note that while the results herein have focused on the 
construction of unitaries, the ideas can be extended to simulat- 
ing non-unitary processes such as generalized measurements. 
Generalized measurements on a system s can be thought of 
as orthogonal measurements on an extended system H s ® H^, 
which may not be orthogonal in s alone. Applications includ- 
ing precision measurement [18], quantum communication in 
the context of entanglement purification [19], and quantum 
error correction [20]. To realize this positive operator val- 
ued measurement (POVM), one can perform a unitary oper- 
ation on H s © Hf- followed by a projective measurement on 
alone. For example, non-orthogonal measurements on a 
qubit can be realized by appending ancillary qubits, perform- 
ing unitary operations on the joint system, and measuring the 
ancillae. The requirement of using two qubit gates can be ob- 
viated if the ancillary degrees of freedom come from orthog- 
onal states within the same system. For example, one can use 
the d — 2 states of a qudit to implement POVMs on an orthog- 
onal qubit subspace. These ideas are explored in the context 
of quantum optical systems in Refs. [21, 22]. The techniques 
reported here indicate that the requisite operations on the ap- 
pended Hilbert space can be done efficiently. 
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APPENDIX A: CONNECTION TO EARLIER MULTIQUDIT 
GATE CONSTRUCTIONS 

Mathukrishnan and Stroud have shown [ ] that exact uni- 
versal computation over qudits can be achieved using a gate li- 
brary containing a 2d — 1 parameter family of two qudit gates. 
We show that this family of gates can be simulated using d — 1 
instances of a single parameter two-qudit gate generated by 
the Hamiltonian H lnt (Eq. 11). 

They begin by writing the unitary W € SU(d") in its spec- 
tral decomposition: 



The unitary can then be expressed as the product 

W = Y\X f (\Xj))V f (Xj)Xf\\Xj)), (A2) 
;'=o 

Here the operator Vf(kj) applies a phase only to a fiducial 
logical basis state |/), 

V f (X J )=e iX J\f){f\+l i \k)(k\, (A3) 

¥/ 

The operator X/( | Xj}) is a unitary extension of the map from 
the fiducial state to an eigenvector of U: 

Xj(\Xj)) = \Xj){f\+Z,\x k (j))(k\. (A4) 

Hi 

where (^|x-t0')> = and (Xk(j)\%k?U)) = 8 W- There is a 
freedom in the choice of the unitary extension by fixing the 
set of mappings {\k} — > \%k(j))}- Notice that arbitrary sin- 
gle qudit operations can be constructed using the spectral de- 
composition for n = 1 and choosing the fiducial state to be 
logical basis state of one qudit. Here we fix |/} = \d — 1). 
The two multiqudit operators Eqs. A3, A4 can be simulated 
exactly the using single qudit operations and two families of 
controlled two-qudit operators. The first two-qudit gate de- 
fines a one parameter family of controlled-phase gates and is 
in fact generated directly by //,„, : 

Ai(V{(^))=e- m ^^ m \ (A5) 

The second family of operators is defined Ai(X(\\\t))) and 
maps \d — 1) — ► |\|/) on the target qudit iff the control is in 
state \d — 1) and applies 1 to the target otherwise: 

Ai(x(| V ») = iw-ij/wm 

+\d-\){d-\\®{\y)(d-\\ (A6) 

+Em*-iIP*X*I)> 

where (y|pt) = and (P/|pi) = S ;J t. Because the gate is al- 
lowed to implement any unitary extension of \y)t(d — 1 1, it 
only depends on the 2d — 2 parameters of the state |\|/) (two 
parameters are fixed by the norm (\|/|\|/) = 1 and setting the 
global phase to zero.) 

The gate A\(X(\\\i))) can be simulated exactly with the 
controlled-phase gate (Eq. A5) and single qudit gates as we 
now show. First, expand the state \\\t) in the single qudit ba- 
sis: |\|/) = TfjZo c j\j)' where the global phase is chosen so that 
argcd-i = 1. The conditional mapping \d — 1} — > can be 
realized as a sequence of d — 1 controlled unitaries that couple 
two target qudit basis states at a time, 

Ai(X(|v))) = nAi(^/-i(y/,^)), (A7) 

.7=0 

The arguments (Qj , <|)y) for each controlled unitary must sat- 
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isfy the following relations: 





2 


Ud-2,d- 


i\d-l) 


(d- 


3 


Ud-3,d- 


[Ud-2.d- 


(d- 


3 


Ud-3,d- 




(d- 


1 


Ud-2,d- 


,\d-l) 


—ie 


-*rf-3 siny f /_ 


-3 COS J d - 



(A8) 



c k = -/e* sin Y*nf=t + i cosy/ (k<d-2). 
Now it only remains to demonstrate that each controlled ro- 



tation Ai(I//d_i) can be simulated with just the controlled 
phase gate and rotations on the target qudit. A single conju- 
gation suffices: 



M{U jjt -i) = KZ>Uj,d-i(yj/2 + K,$j) 

(A9) 

Following this construction, d — 1 controlled phase gates and 
d single qudit gates suffice to exactly simulate Ai )). 
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